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ABSTRACT 

In this paper we study the evolution of ^Li in the Galaxy considering the con- 
tributions of various stellar sources: type II supernovae, novae, red giant stars, 
and asymptotic giant branch (AGB) stars. We present new results for the pro- 
duction of ''Li in AGB stars via the hot bottom burning process, based on stellar 
evolutionary models by Frost (1997). In the light of recent observations of dense 
circumstellar shells around evolved stars in the Galaxy and in the Magellanic 
Clouds, we also consider the impact of a very high mass-loss rate episode (su- 
perwind) before the evolution off the AGB phase on the ^Li enrichment in the 
interstellar medium. We compare the Galactic evolution of ^Li obtained with 
these new ''Li yields (complemented with a critical re-analysis of the role of su- 
pernovae, novae and giant stars) with a selected compilation of spectroscopic 
observations including halo and disk field stars as well as young stellar clusters. 
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We conclude that even allowing for the large uncertainties in the theoretical cal- 
culation of mass-loss rates at the end of the AGB phase, the superwind phase 
has a significant effect on the ^Li enrichment of the Galaxy. 

Subject headings: nucleosynthesis - stars: abundances, AGB and post-AGB - 
Galaxy: evolution, abundances 

1. Introduction 

In spite of many attempts, the evolution of ^Li in the Galaxy is not completely un- 
derstood. Lithium has the unique property of being produced by, at least, three different 
processes: Big Bang nucleosynthesis, spallation of Galactic cosmic-ray (GCR) particles on in- 
terstellar matter (ISM) nuclei, and stellar nucleosynthesis. At the same time, ^Li is easily de- 
stroyed by proton captures in stellar interiors at relatively low temperatures (~ 2.5 x 10^ K). 

Observationally, Population I (hereafter Pop. I) dwarf stars show a large dispersion in 
^Li abundance attributed to different amounts of depletion. On the contrary the abundance 
of ^Li in the atmosphere of old, warm (Tcfj > 5700 K) Population II (hereafter Pop. II) 
dwarf stars is surprisingly uniform: ^Li/H ^ 1.1 x 10~^° (Spite & Spite 1982), a value 
~ 10 times lower than the maximum abundance observed in Pop. I stars. Several studies 
have subsequently confirmed the existence of the so-called "Spite plateau" in Pop. II stars, 
although discrepancies on the actual value of the plateau still exist at the level of ~ 0.1 dex 
(see e.g.. Spite & Spite 1986; Hobbs & Duncan 1987; Rebolo, Beckman, & Molaro 1988; 
Thorburn 1994; Spite et al. 1996; Bonifacio & Molaro 1997) Q 

There is an ongoing debate about whether the Spite plateau corresponds to the (al- 
most) undepleted primordial abundance of ''Li , or whether the primordial value is closer to 
loge(^Li) ~ 3.1-3.3 (the value measured in the youngest stars and solar system meteorites) 
and has been significantly but uniformly depleted in Pop. II stars (see e.g. Deliyannis & Ryan 
1997; Bonifacio & Molaro 1997). On one hand, the former interpretation requires processes 
that increase the Galactic ''Li abundance by a factor ~ 10 in the first few Gyr. On the other 
hand standard models of stellar evolution which incorporate only convective mixing predict 
little (~ 0.1 dex) ^Li depletion in Pop. II stars (e.g., Deliyannis, Demarque & Kawaler 1990). 
Stellar models including the effects of non-standard processes (like e.g. diffusion, winds, tur- 
bulence induced by rotational instabilities, slow mixing driven by angular momentum loss. 



^In this paper we adopt the usual notation loge(^Li) — log[7V(''Li)/A^(H)] + 12. The Spite plateau 
corresponds to loge(^Li) ~ 2.1. 
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etc.) predict substantial ^Li destruction in Pop. II stars (see Pinsonneault 1997 for a review). 
These non-standard models predict observable features, like a dispersion in ''Li among halo 
dwarfs, or correlations between ^Li and Teg or [Fe/H], which have not been confirmed by 
observations, suggesting that ^Li destruction must have been minimal in these stars (see 
e.g., Boesgaard & Steigman 1985, Bonifacio & Molaro 1997, Ryan, Norris, & Beers 1999, 
Ryan et al. 2000). It is true, however, that the presence of a few metal-poor stars with 
^Li abundances somewhat above the plateau, as well as other observational evidences in 
Pop. I stars (Deliyannis et al. 1998, Deliyannis 2000) support the idea that halo dwarfs may 
have suffered a certain amount of ^Li depletion. 

Indication of a primordial abundance of ^Li in excess of the Spite plateau value is sug- 
gested also by (i) recent determinations of deuterium abundances in Lyman-a clouds (see 
e.g. Pettini & Bowen 2001 and references therein), and {ii) observations of the spectrum of 
fluctuations in the cosmic background radiation (de Bernardis et. al. 2000, Tegmark & Zal- 
darriaga 2000). In both cases the value of the baryon-to-photon ratio r/io results in the range 
6-7, rather than 4-5 as usually assumed. This implies a cosmological abundance of ^Li in 
excess of the Spite plateau value by a factor 2-3, still compatible with the observationally 
allowed range of primordial ^He. However, one should be very careful in interpreting these 
indications since possible sources of systematic errors have still to be completely understood 
and eliminated from the observational data. The issue remains therefore open. 

Assuming that the Spite plateau corresponds to the primordial ''Li abundance, the 
difference between the plateau abundance and the maximum ^Li content in Pop. I stars 
underscores the need for one or more sites of ^Li production. 

It is well known that ^Be, ^"-^^B, and ^-'^Li can be produced via spallation and fusion 
reactions between GCR particles and ISM nuclei (Fowler, Reeves, & Silk 1970; Meneguzzi, 
Audouze & Reeves 1971). Whereas theoretical calculations of the amounts of ^Be, ^^'^^B and 
^Li produced by GCR (see e.g. Lemoine, Vangioni-Flam, & Casse 1998; Fields & Olive 1999; 
Ramaty et al. 2000, and references therein) can roughly reproduce the observed values, the 
amount of ^Li produced by spallation processes at the time of formation of the Sun is a factor 
5-10 lower than the measured meteoritic value (see also Nissen et al. 1999 and references 
therein)^. 

Additional sites of ^Li production are therefore required to match the observational 



^At lower metallicity, in particular in the range — 2 < [Fe/H] < —0.5, the predictions of GCR spallation 
models indicate a possible overproduction of ^'^Li due to a-a reactions (Prantzos et al. 1993, Valle et al. 2001 
in preparation) that may suggest the occurence of a depletion process in halo stars at least in this metallicity 
range. 
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data. For example, D'Antona & Matteucci (1991) proposed that novae and AGB stars may 
represent the main sources of Galactic ^Li , whereas Matteucci, D'Antona, & Timmes (1995) 
favoured a model where the main contributors to ^Li enrichment are SNll and AGB stars. 
Abia, Isern, & Canal (1995) considered the production of ^Li by low-mass AGB stars, whereas 
Romano et al. (1999) showed that novae are required in order to reproduce the growth of 
^Li with metallicity, and that the contribution from SNII should be lowered by at least 
a factor of two. Finally, Romano et al. (2001), using recent AGB models by Ventura, 
D'Antona, & MazziteUi (2000), concluded that AGB stars are not substantial ^Li producers 
in the Galaxy. 

In this paper we reanalyze the problem afresh with the help of a numerical model of 
Galactic chemical evolution and a set of new models of AGB stars of various masses and 
metallicities. We consider various stellar sources of ^Li : novae, type II supernovae, low- 
and intermediate-mass stars in the red giant branch (RGB) and AGB phase. All these 
sources can produce ^Li , as shown by theoretical models, and, in some cases, also confirmed 
by spectroscopic data. However, individual ^Li yields are rather uncertain and strongly 
dependent on model parameters and assumptions. 

The paper is organized as follows. In § 2 we discuss our sample of observational data for 
field stars and open clusters. In § 3 we briefly describe the adopted model for Galactic chem- 
ical evolution. In the following Sections we then focus on the analysis of the ^Li production 
by different stellar sources. In particular in § 4 we consider the role of intermediate-mass 
AGB stars (hereafter IMS-AGB), presenting new ^Li yields obtained from the AGB models 
by Frost (1997). We reanalyze the production of ''Li by novae (§ 5), SNII (§ 6), and low-mass 
RGB stars (§ 7), on the basis of recent model calculations and observations. Finally, in § 8 
we discuss and summarize our results. 



2. Observational data 

2.1. ^Li abundances in field dwarfs 

The main source of our compilation is the recent study by Fulbright (2000) who de- 
termined ^Li , Fe, and other element abundances for 168 halo and disk stars. Both stellar 
parameters and abundances of ^Li and Fe were inferred by carrying out a self-consistent LTE 
analysis, providing a homogeneous sample of both Pop. II and Pop. I stars. The survey of 
Fulbright does not extend to very low or very high metallicities; therefore we complemented 
it for [Fe/H] < -2 with the compilation of Bonifacio & Molaro (1997). At [Fe/H] > -0.4 
we included a sample of warm (Teff > 5700 K), unevolved F-type stars taken from Balachan- 
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dran (1990) and Lambert, Health, & Edvardsson (1991). We note that, whereas cooler stars 
have certainly undergone Li depletion (unless very young), considering only stars warmer 
than 5700 K allows having a subsample of minimally depleted stars which trace the upper 
envelope of the ^Li vs. [Fe/H] distribution. Surveys of ^Li among Pop. II stars have also been 
made by e.g. Rebolo et al. (1988), Thorburn (1994), Ryan, Norris, & Beers (1996), Spite et 
al. (1996), Ryan et al. (2001). Different datasets and/or analysis methods (in particular the 
use of different effective temperature calibrations) lead in general to a different abundance 
scale: therefore, including ^Li abundances from many different sources in the same compila- 
tion may introduce a spurious (i.e. not real) dispersion in the ^Li vs [Fe/H] distribution. For 
this reason, we decided to base our compilation of ''Li abundances for Pop. II stars on the 
minimum number of studies allowing the full coverage of the entire metallicity range. We 
mention that Fulbright (2000) derived an average offset between his and Bonifacio & Mo- 
laro's (1997) effective temperatures ATgs = —38 ± 20 K, which should not introduce major 
systematic differences in the inferred ^Li abundances. There are several stars in common 
between the two studies: for all of them the agreement in ^Li abundances is indeed very 
good (within ~ 0.1 dex). 

We did not include in our list stars evolved off the main sequence (hereafter MS) and/or 
stars cooler than 5700 K; both groups of stars most likely did not preserve their original 
^Li content due to either post-MS dilution or to MS ^Li depletion and thus are not adequate 
to trace the upper envelope of the ''Li vs. metallicity distribution. We also did not include 
''Li abundances for T Tauri stars, since they may be affected by several problems such as 
spectral veiling, NLTE, or, more generally, by the effects of a "disturbed stellar photosphere" 
(see e.g. Duncan 1991). Instead, we included in our compilation the average abundance of 
warm (i.e. F-type) stars in young open clusters with —0.1 <[Fe/H]< 0.1. These stars have 
undergone minimal pre-MS or MS ^Li depletion and thus their abundance can be considered 
as representative of the present ISM ^Li content. Several surveys of ^Li among open cluster 
stars (see e.g. Jeffries 2000 for a recent review) have shown that Li-undepleted stars in all 
the investigated clusters share the same abundance loge(''Li) = 3.2 ± 0.1 (e.g., Randich et 
al. 1997)0. 

In the upper panel of Figure 1 we plot the ^Li abundance vs. metallicity for the stars 
of our sample. A few features in this figure warrant detailed comments. As far as the 
''Li plateau is concerned, recent studies have inferred slightly different values for the plateau 



■^Although ^Li abundances as high as loge(^Li) = 4 have been derived for some T Tauri stars, Magazzu, 
Rebolo, & Pavlenko (1992) and Martin et al. (1994) concluded that abundances of pre-MS stars (including 
both classical and weak-lines T Tauri stars) as a whole indicate an initial ''Li abundance for Pop. I stars 
loge(^Li) = 3.1. 
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value, due - as mentioned above - to different types of analysis and, in particular, to possible 
offsets in the zero-point of the effective temperature scale (see e.g., Spite et al. 1996). In our 
model of Galactic evolution we have generally adopted as initial ^Li abundance the value 
determined by Bonifacio & Molaro (1997) for the Spite plateau, i.e. loge('^Li)o = 2.24±0.012, 
under the implicit assumption that the Spite plateau represents the primordial (Big-Bang) 
abundance. We mention in passing that the mean abundance of the stars in the compilation 
of Fulbright (2000) with [Fe/H]< -1 is loge(^Li)o = 2.24 ± 0.095 in excellent agreement 
with the estimate of Bonifacio & Molaro (1997). We also ran our model starting with an 
initial ^Li abundance higher by a factor 2 than the Spite plateau value (see below) , following 
the idea that Pop. II stars may have depleted a certain amount of Li in the course of their 
evolution. 

Figure 1 also suggests that the rise of ^Li abundance from the Spite plateau may be 
smoother than claimed in previous studies. In particular, a few stars are present in the survey 
of Fulbright (2000) which have intermediate metallicities (—0.7 < [Fe/H] < —0.45) and 
''Li abundances 0.2-0.3 dex (i.e. > la) above the plateau but not as high as loge(^Li) ~ 3. 
However, the caveat is that only a few stars define the upper envelope at intermediate 
metallicity and we cannot exclude that these stars have suffered some Li depletion. 

A sharp increase to abundances as high as the present ISM abundance may occur at 
metallicities between [Fe/H]= —0.45 and —0.35, if the abundances for the few stars which 
trace the upper envelope of the distribution at those metallicities are correct. These stars 
are from the survey of Balachandran (1990), but we mention that Li abundances derived 
by Balachandran are systematically higher than those of Boesgaard & Tripicco (1986), with 
differences in the range 0.2 - 0.6 dex. For the case of HR 8315 (HD 206901) Lebre et 
al. (1999) inferred a metaUicity [Fe/H]= —0.3 and a ^Li abundance loge(Li) = 2.9 to be 
compared with [Fe/H]= —0.37 and loge(Li) = 3.05 determined by Balachandran (1990). We 
also notice that photometric iron abundances for some of the intermediate metallicity and Li- 
rich stars observed by Balachandran (1990) are at least 0.2 dex higher than the spectroscopic 
values inferred by Balachandran (see the last two columns of Table 1 in Balachandran 1990), 
suggesting again that ^Li abundance ~ 3 is reached only at higher metallicities (i.e. [Fc/H]> 
—0.3). Determinations of ''Li abundances in this critical metallicity range are certainly 
needed in order to better constrain the ''Li vs. [Fe/H] upper envelope. 

2.2. ^Li abundances in clusters 

In the lower panel of Figure 1 we show the abundance of ''Li vs. Galactic age us- 
ing "^Li data for stars in Galactic open and globular clusters. We also show the meteoritic 
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^Li abundance (loge(^Li) = 3.31 ± 0.04, Anders & Grevesse 1989), the interstellar ^Li abun- 
dances determined in the line-of-sight toward p Oph (^Li /H ^ 2.4 x 10~^°, Lemoine et 
al. 1993), and toward o and C Per f Li /H = (9.8 ± 3.5) x IQ-^o and (12.2 ± 2.2) x 10"^°, 
Knauth et al. 2000). For young open clusters (ages < 1 Gyr) the average abundance mea- 
sured in warm (Tcfj > 6000 K), undepleted cluster MS stars, is plotted in the Figure. The 
data were taken from the following sources: Sodcrblom et al. (1999) for NGC 2264; Randich 
et al. (1997) for IG 2602; Randich et al. (2001) for IG 2391; Martm & Montes (1997) for 
IG 4665; Randich et al. (1998) for Alpha Persei; Soderblom et al. (1993) for the Pleiades; 
Thorburn et al. (1993) for the Hyades. On the other hand, the initial abundance for older 
clusters was estimated from stars at (or just evolved off) the turn-off in the 2 Gyr clusters 
NGG 752 (Hobbs & Pilachowski 1986), NGG 3680 (Randich, Pasquini, & Pallavicini 2000) 
and from the abundances of the two components of the tidally locked binary S1045 in the 
solar-aged M 67 cluster (Deliyannis ct al. 1994; Pasquini, Randich, & Pallavicini 1997). Fi- 
nally, globular clusters (GCs) arc represented by a box that shows the ranges in '^Li and age 
estimates. More specifically, GCs absolute ages are the subject of a long standing debate; in 
the lower panel of Figure 1 we used the recent determination by Garretta et al. (2000), who 
inferred an average age of 12.9 ± 2.9 Gyr based on Hipparcos parallaxes for local subdwarfs. 
We mention, however, that ages as large as 16-17 Gyr cannot be ruled out. 

Abundances of ^Li were derived for unevolved stars at the turn-off in three clusters only, 
namely, NGG 6397 (Pasquini & Molaro 1996), 47 Tuc (Pasquini & Molaro 1997), and M 92 
(Boesgaard et al. 1998); whereas the average abundances for the three clusters are close 
to the value of the Spite plateau, stars in the same cluster show a significant dispersion in 
^Li (see discussion in Boesgaard et al. 1998) and, in particular, a star is found in M 92 with 
a ^Li content significantly higher than the plateau and the other chistcr stars. Consequently, 
we think that the use of a box is more appropriate to indicate the position of GCs in Figure 1. 



3. The Galactic chemical evolution model 

The model of Galactic chemical evolution adopted in this work is described in detail by 
Ferrini & Galli (1988), Galli & Ferrini (1989) and Ferrini et al. (1992). The same model was 
adopted by Galli et al. (1995), Travaglio et al. (1999) and Travaglio et al. (2001), to study 
the evohition of the light elements D and '^Hc, and the evolution of the heavy elements from 
Ba to Pb, respectively. We briefly recall here the basic features of the model. 

The Galaxy is divided into three zones, halo, thick disk, and thin disk, whose composi- 
tion of stars, gas (atomic and molecular) and stellar remnants is computed as functions of 
time up to the present epoch tgai = 13 Gyr. The thin disk is formed from material infalling 
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from the thick disk and the halo. The formation of the Sun in the thin disk takes place 
4.5 Gyr ago, i.e. at epoch = 8.5 Gyr. The star formation rate in the three zones is not 
assumed a priori, but is obtained as the result of self-regulating processes occurring in the 
molecular gas phase, either spontaneously or stimulated by the presence of other stars. 

In models of Galactic chemical evolution it is customary to separate the contribution of 
stars in different mass ranges to the enrichment of the ISM. Our model follows the evolution 
of (i) single low- and intermediate-mass stars (0.8 Mq < M < M^) ending their life as He 
or C-0 white dwarf, (ii) binary systems able to produce type I supernovae, and {Hi) single 
massive stars (M^ < M < 100 Mq), the progenitors of type II supernovae. The value of 
depends on metallicity; we assume = 6 Mq for [Fe/H]< —0.8, and M^ = 8 Mq otherwise 
(see for references Tornambe & Chieffi 1986). The adopted initial mass function (hereafter 
IMF) is discussed in Ferrini et al. (1992). Stellar nucleosynthesis is treated according to 
the matrix formahsm introduced by Talbot & Arnett (1973). Specific yields are taken from 
calculations by Woosley & Weaver (1995) and Thielemann, Nomoto, & Hashimoto (1996), 
for type II and type I supernovae, respectively. 

4. The production of '^Li by hot bottom burning process 

Hot bottom burning (hereafter HBB) occurs in intermediate mass AGB stars when 
the bottom of the convective envelope of the star reaches the top layers of the hydrogen 
burning shell. Sackmann & Boothroyd (1992) showed that the Cameron-Fowler mechanism 
(Cameron 1955; Cameron & Fowler 1971) is very effective at producing ^Li , and can quanti- 
tatively account for the ^Li abundance measured by Smith & Lambert (1990) in the super-Li 
rich AGB stars in the Magellanic Clouds. We present here a summary of results for ^Li pro- 
duction in AGB stars of various masses and metallicities, based on the evolutionary models 
of Frost (1997) and synthetic models using the technique of Forestini & Charbonnel (1997). 
We then discuss the sensitivity of ^Li enrichment to the mass-loss prescriptions during the 
latest phases of AGB evolution. 

4.1. Li production from 4, 5, 6 Mq AGB models at different metallicities 

Frost (1997) computed the evolution of 4, 5, and 6 Mq models each with metal content 
Z — 0.02, 0.008, and 0.004, appropriate to the solar neighborhood, and to the Large and 
Small Magellanic Clouds, spanning the evolution from the MS through the end of the AGB 
phase. Each case has been the subject of a detailed nucleosynthesis study using a modified 
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post-processing code (Cannon 1993) including 74 species and over 500 reactions. Some of 
the results of these models have been already published (e.g. Frost et al. 1998, Lattanzio & 
Forestini 1999); in this paper we present in detail an analysis of the production of '^Li. 

Our present approach is to use the stellar models of Frost (1997) as the basis of synthetic 
models computed with the code of Forestini & Charbonnel (1997). In this way we force the 
synthetic code to follow the results of the stellar models (accuracy is within 10%), yet we 
are free to vary those parameters which have no significant feedback on the evolution. In 
particular, for the aim of this paper, we investigated the mass-loss rate at the end of the 
AGB phase (see next Section for detailed discussion). We do not expect modest changes in 
the mass-loss rate during the latest AGB phase to have a significant effect on the structural 
evolution of the star, since the variations in the mass-loss rate only affect the very late 
phases of stellar evolution. Nevertheless, we caution against placing too much behef in the 
quantitative results for the extreme cases discussed below. 

The surface ^Li abundances during the AGB lifetime in the stellar envelope for the 
nine models considered here is presented in Figure 2. Note that the abundance of ^Li at 
the beginning of the AGB phase varies from one model to another as a result of previous 
evolution. Fig. 2 is a 3 x 3 grid figure with mass increasing along the x-axis and Z decreasing 
along the y-axis: in this way we expect HBB to increase for increasing x-axis (mass) or 
increasing |/-axis (decreasing Z\ The most extreme HBB is for the top right plot, i.e. 
M = 6 Mq model with Z = 0.004, where the maximum ^Li production lasts for the longest 
time (~ 10^ yr). The maximum '^Li abundances we obtain range around logeC^Li) ~ 4, 
independent of stellar mass and metallicity, and only depending on the occurrence of HBB. 
The models with M = 4 Mq and Z = 0.02, 0.008 show no HBB and consequently no 
^Li production. In fact, ^Li destruction in these two models proves that their convective 
bottom temperatures are high enough to burn ^Li through ^Li(p, ci;)^He channel, but not 
enough to produce ^Be (that decays to ^Li, Cameron-Fowler mechanism) by ^He burning. 
Conversely the 6 Mq models, with their substantial envelopes, begin HBB immediately on 
the AGB and hence they produce a large amount of ^Li , but the high envelope temperature 
also assists in the destruction of ^Li . Finally, Figure 2 shows that in all models, as the 
evolution of the star proceeds, the surface ^Li abundance decreases again more or less rapidly, 
depending on the mass and metallicity. This indicates that the initial reserves of '^He have 
been consumed in the production of ^Li , and this ^Li is now itself being destroyed (without 
replenishment) by proton captures. 

There are observational indications that some Galactic C-stars (stars with C/0 > 1) 
are also super-Li rich stars (loge('^Li) > 3, see e.g. Abia et al. 1991, Abia, Pavlenko, & de 
Laverny 1999, Abia & Isern 2000). Therefore, we also examined the ^Li vs. C/0 trends 
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for the nine cases presented above. The models of Frost (1997) use the algorithm for the 
third dredge-up (i.e. the penetration of the convcctivc envelope into the partially He-burnt 
zone after each thermal pulse), as described by Frost & Lattanzio (1996), and include the 
entropy adjustment by Wood (1981). After each third dredge-up episode the ^^C abundance 
is increased. Hence third dredge-up increases the ratio C/0 (by adding C to the stellar 
envelope), and HBB (where it is active) decreases it again as it transforms C into N. Figure 3 
shows the C/0 ratio vs. ^Li for the nine models under consideration. In some cases HBB 
prevents the formation of a C star, while in other cases it merely delays the formation to 
a time when HBB has stopped and third dredge-up continues (see Frost et al. 1998 for 
more details). The condition C/0 > 1 together with high ^Li values are obtained only 
for the case Z = 0.004 and M = 4 Mq and for a very short phase for the case with the 
same metallicities and M = 5 Mq. In particular, the 4 Mq model is a Li-rich star (with 
loge(''Li) ~ 4) and C-rich star for a period of about 10^ yr. Moreover, all the other cases 
show high ^Li abundances together with 0-rich envelopes, in agreement with observations 
of AGB stars in the Magellanic Clouds (see e.g. Smith et al. 1995) and in the Galaxy (see 
e.g. the recent work from Arellano Ferro, Giridhar, & Mathias 2001). 

The range of Mboi covered by the nine models presented here (—7 < Mboi < —6) during 
the Li-rich phase is in agreement with bolometric magnitudes of the super-Li rich stars 
observed in the Magellanic Clouds (Smith et al. 1995). Concerning Galactic super- Li rich 
C stars, only few estimates of their Mboi are available. Abia et al. (1991), based on the 
distance determinations by Claussen et al. (1987), estimated Mboi — —5 for WZCas and 
WXCyg (with loge(^Li) = 5.0 and 4.7, respectively). More recently, Abia & Isern (2000), 
based on Hipparcos parallaxes, inferred Mbd = —6.44 and loge('''Li) =4.8 for WZCas, and 
Mboi = —4.35 and logeC^Li) = 4.4 for WXCyg. However, a caveat is that the latter is 
based on a parallax vr = —1.41 ± 1.98, which is unreliable. Two additional Galactic C-rich 
super-Li rich stars have been presented by Abia et al. (1991): lYHya (Mboi = —6.2 and 
loge(^Li) = 5.4), and TSgr (Mboi = —5.8 and loge(^Li) = 4.2). Given the uncertainties, 
and excluding the not rehable measurement of Mboi for WZCas by Abia & Isern (2000), we 
conclude that the Mboi determinations for the Galactic super-Li rich stars are still consistent 
with the Mboi reached by our AGB models during the super-Li rich phase. We also notice 
that the highest abundance reached by our models for log e(^Li) is ~ 4, while higher values 
are inferred from observations. 
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4.2. Superwind in AGB stars: consequences for Li chemical enrichment 

There are observational indications that toward the end of their evolution, AGB stars 
(both of low- and intermediate-mass) can experience a short phase of extremely rapid mass- 
loss, called superwind (van der Veen, Habing, & Geballe 1989; Schroder, Winters, & Sedl- 
mayr 1999). In this phase periods of heavy mass-loss may be interspersed with much longer 
periods of lower mass-loss rates, possibly powered by He shell flashes. These stars are of- 
ten surrounded by a dense circumstellar envelope and some of them are no longer visible 
at optical wavelengths. Evidence for a superwind has been found in CO observations of C 
stars (e.g. Olofsson 1993), of K giants (de La Reza et al. 1997; Castilho et al. 1998), and in 
post-AGB objects that show concentric rings around the central object (see e.g. Crabtree & 
Rogers 1993; Latter et al. 1993; Riera et al. 1995; Klochkova et al. 1999; Blocker et al. 2001). 

From the theoretical point of view, the term of superwind has been coined (Renzini 1981) 
to describe the heavy, final tip- AGB mass-loss (> 10"^ Mq yr^^), which is required to form 
a planetary nebula of typically tenths of a solar mass within several 10^ yr. The idea of 
a star terminating its AGB life with a superwind has been also elaborated by Bowen & 
Willson (1991). These authors showed that all stars undergo the superwind phase as a result 
of an increase in scale height and density at the condensation radius as the star evolves 
toward the tip of the AGB. Analytical expressions to describe this phase of high mass- 
loss have been presented by different authors, e.g. by Vassihadis & Wood (1993, hereafter 
VW93), Blocker (1995), Salasnich, Bressan, & Chiosi (1999). VW93 computed the effects of 
thermal pulses on mass-loss rates and suggested that a star may undergo several superwind 
phases and that more massive stars will stay longer in the superwind phase than lower mass 
stars. For the AGB models discussed in this Section the superwind prescription follows the 
formula of VW93 (without the correction for masses above 2.5 Mq) 



where the mass-loss rate is in Mq yr ^ and P is the pulsation period in days. Note that the 
mass-loss rate is truncated at 



where L is the stellar luminosity, corresponding to a radiation-pressure driven wind. Again 
from VW93 we take the wind expansion velocity (in km s~^) 




(1) 



^;exp = -13.5 + 0.056P, 



(3) 



ranging between 3.0 and 15.0 km s 



-1 
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The ^Li yields are very sensitive not only on the extent of the HBB, but also on the 
mass-loss rate prescriptions. As first qualitatively pointed by Abia, Isen, & Canal (1993), 
it is not clear if during the super-Li rich phase AGB stars are able to inject ''Li into the 
interstellar medium by the high mass-loss rate before ''Li is depleted in the atmosphere. 
We analyze this point with our models, and we plot in Figure 4 the surface ''Li abundance 
as a function of the mass-loss rate for the nine models described in the previous Section. 
The figure shows that only in a few cases does the star achieve very high mass-loss rates 
(> 10~^ Mq yr~^) when the surface ^Li abundances is still high (loge(^Li) > 3). When this 
happens the overall ^Li yield is positive, or only slightly negative. In the first three columns 
of Table 1 we report the ^Li yields (in Mq) obtained with these models, showing that only 
in the case of 4 Mq and Z — 0.004 the ^Li yield is positive. 

Another way to investigate the sensitivity of ^Li yields to the choice of the mass-loss rate 
is shown in Figure 5, where we plot the surface ^Li abundance as a function of the total stellar 
mass (so that the models evolve from right to left). Note that for the models with Z — 0.02 
we assumed an initial (i.e. at the time when the star formed) loge(^Li) ~ 3.3, compatible 
with the meteoritic abundance; this value scales with Z for the other metallicities. In order to 
get a positive yield of ^Li the star's mass should decrease substantially while loge(^Li) > 3-4. 
Hence, although the 5 and 6 Mq models experience HBB, the 5 Mq loses about 1 Mq when 
log e(^Li) ~ 4, whereas the 6 Mq loses most of its mass when log e(^Li) is 3.5 or lower. Hence 
the yield from 5 Mq model, although still negative, is smaller in its absolute value than the 
6 Mq model. That is, the 5 Mq model does not destroy as much '''Li as the 6 Mq model, 
due to the large amount it returns to the ISM, although the overall yield is still negative, 
with respect to the initial abundance. At a lower value of Z (~ 0.008) where typical interior 
temperatures are higher, both the 5 Mq and 6 Mq models clearly show HBB, and there is a 
small amount of HBB at 4 Mq, although loge(^Li) never exceeds 2. Again, most of the mass- 



TABLE 1 

'''Ll YIELDS (in solar MASSES) FROM INTERMEDIATE-MASS AGB STARS 





(a) 


(b) 




Z = 0.004 


Z = 0.008 


Z = 0.02 


Z = 0.004 


Z = 0.008 


Z = 0.02 


4M0 


1.1 X 10-^ 


-1.7 X 10-^ 


-4.3 X 10-^ 


1.5 X 10-^ 


-1.3 X 10-s 


-3.3 X 10-s 


5 Mq 


-1.1 X 10"^ 


-1.9 X 10-s 


-4.0 X 10-9 


4.5 X 10-*^ 


7.3 X 10-s 


2.2 X 10-8 


6 Mq 


-1.3 X 10"^ 


-2.5 X 10"^ 


-1.3 X 10"^ 


4.1 X 10"^ 


7.2 X 10"^ 


7.7 X 10-8 



(a) Mass-loss rate from VW93 

(&) Mass-loss rate from VW93 increased by a factor 50 
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loss for these models occurs when the surface e(^Li) is less than the initial value, and hence 
the overall yields arc negative (see Table 1). Finally, at Z = 0.004 the higher temperature 
means that there is significant ^Li production even at 4 Mq. This is particularly strong in 
the more massive models, but this does not correspond to a positive ^Li yield because the 
^Li is destroyed before the high mass-loss begins. In contrast, the 4 Mq model has its HBB 
delayed, and toward the end of its evolution the increased mass-loss begins when the surface 
^Li is higher, and hence the overall yield is higher. Also note that the high values found 
for log e(^Li) in the 5 and 6 Mq models do not last until the mass-loss reaches appreciable 
values. By the time that the mass-loss begins in earnest, the surface ^Li abundance is below 
the initial value. For the 4 Mq case, however, there is about a half solar mass of material 
lost when loge(^Li) is above 3. A conclusion from these results is that the higher ^Li yields 
come from the lower mass stars: in fact, from those stars which just start HBB when the 
high mass-loss rate begins. 

Due to the fact that ^Li yields are very sensitive to the mass-loss choices, we have run 
each of the models presented above using different mass-loss rates. Since several observations 
of OH/IR stars (i.e. 0-rich AGB stars that exhibit OH masers, Wilson & Barrett 1972) with 
infrared excesses (i.e. high mass-loss rates, see e.g. Blocker et al. 2001) and P ~ 400-500 
days are now available (see e.g. Lewis 2000 for a recent survey of OH/IR IRAS sources), we 
forced the mass-loss to start at shorter periods. In order to obtain dM/dt > 10~^ Mq yr~^ 
with P ~ 500 days the VW93 prescription has to be increased of a factor of ~ 50. In Table 1 
(columns 5, 6, 7) we hst the yields obtained with this modification to the VW93 prescription. 
The surface ^Li abundance for the different models vs. stellar mass is also shown in Figure 5 
for comparison. Both the table and the figure indicate that, when using this mass-loss rate 
prescription, the ^Li yields for most of the cases are positive. 

The chemical evolution model described in § 3 has been run using the Li yields given 
in Table 1 (first three columns): we obtain that intermediate- mass AGB stars contribute 
to the solar system ^Li abundance for ~ 14%. Under these conditions, IMS-AGB stars can 
contribute only a small fraction to the ''Li solar composition and they do not seem able to 
reproduce the rapid increase of ''Li in the Galactic disk. When we introduce in the GCE 
model the ''Li yields obtained the modified VW93 mass-loss prescription (Table 1, last three 
columns) and we find that the contribution of intermediate-mass AGB stars to the solar 
^Li composition increases by a few percent (~ 6%). In Figure 6 {upper panel) we compare 
the observed loge(^Li) vs [Fe/H] distribution with the predictions of the model including 
the modified VW93 mass-loss prescription. 

We also notice that we hmited the IMS-AGB upper mass to 6 Mq . This is due to the 
fact that we did not yet run models for 7 and 8 Mq AGB stars.. Nevertheless we extrapolate 
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the 6 Mq yields to 7 and 8 Mq models to test the 4-8 Mq IMS-AGB mass range. We found 
that the results change only by few percent, due to the lower weight on the IMF of the 7 
and 8 Mq with respect to the 4 Mq stars. 

As a test on the results presented in this Section, we also run several AGB models with 
a Reimers mass-loss (Reimers 1975) and different values of the Reimers parameter rj {rj — 1, 
5, 10 for all masses and metallicities). We find that, even for the higher rj values, the super-Li 
rich phase is reached with a mass-loss rate smaller than ~ 5 x 10^^ Mq yr~^. Therefore, the 
resulting '^Li yields arc very close to those obtained with the VW93 standard case (shown in 
the first three columns of Table 1). However, we stress that large variations in the parameter 
Tj may induce significant feedback on the AGB evolution, and the results of our synthetic 
models may not be completely reliable. 



4.3. Galactic Li enrichment: HBB in 3 Mq AGB stars with low Z? 

It is commonly believed that HBB only occurs for masses greater than about 4 Mq. 
This is not strictly true: the temperature at the base of the envelope has also a strong 
dependence on the metallicity of the star. When the temperature at the base of the deep 
convective envelope is larger than ~ 2 x 10^ K, ^Be is efficiently produced and the surface 
^Li abundance increases. For the models presented here, HBB occurs with initial masses 
> 5 Mq at each metaUicities and > 4 Mq for Z — 0.004. In principle, at lower metaUicities 
the inner envelope becomes hot enough to start HBB at lower masses, but models of lower 
metaUicities through the AGB phase are rare. 

Preliminary calculations for low values of Z have been presented by Lattanzio et al. (2001). 
They found that, for Z — 0.0001, envelope temperatures as high as ~ 2 x 10^ K are reached 

for masses as low as 2.5 Mq. Because of the shape of the IMF we would thus expect these 
stars to contribute substantially to ^Li production, via HBB. In addition, as shown above, 
the maximum abundances obtained in the models discussed in this paper are loge(^Li) ~ 4, 
independently on the stellar mass and metallicity. Therefore, under these preliminary indi- 
cations, we tested the sensitivity of ^Li enrichment when we extend the HBB mass range to 
3-6 Mq at low metallicities {Z < 0.004). For M < 4 Mq we adopted the same yields of the 
4 Mq model presented above (detailed models are in preparation), and for M > 4 Mq we 
used the yields obtained with VW93 increased by a factor of 50 (see Table 1 and discussion 
in the previous Section. The result for GCE of ^Li is shown in the lower panel of Figure 6. 
There is a significant increase of the ^Li enrichment mostly due to the higher weight of these 
stars on the IMF, as well as an important increase of the ^Li contribution of these stars at 
the epoch of the solar system formation (~ 40%). 
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5. The production of "Li by novae 

Novae can in principle contribute to the ^Li enrichment of the ISM via the production 
and subsequent decay of ^Be (i.e. Cameron- Fowler Beryllium mechanism). This process 
can be investigated in detail only with the help of hydrodynamical simulations of the ac- 
cretion and explosive phases of evolution of a nova system. Early works by Arnould & 
N0rgaard (1975), Starrfield et al. (1978), Shara (1980), and Iben & Tutukov (1984) showed 
that ^Li can be produced in considerable amounts during nova outbursts. Subsequent stud- 
ies by Boffin et al. (1993) and Coc et al. (1995) pointed out the sensitivity of the results to 
the adopted nuclear reaction network and the treatment of convection between the accreted 
envelope and the underlying white dwarf core. A systematic analysis of the chemical com- 
position of the ejecta of both CO and ONe novae has been recently presented by Hernanz et 
al. (1996) and Jose & Hernanz (1998, hereafter JH98) by means of hydrodynamical simula- 
tions following both the accretion and the explosion phase. Their results can be summarized 
as follows: (i) the production of ^Be is weakly dependent (logarithmically) on the initial ^He 
concentration in the envelope, for fractional abundances of ^He larger than the solar value; 
(ii) the final ^Li abundance depends rather sensitively on the chemical composition of the 
envelope, which, in turn, depends on the composition of the underlying core: typically, the 
^Li abundance by mass is in the range ~ 10^^-10^^ in the case of a CO white dwarf, and 
~ 10^^-10^® in the case of ONe white dwarfs; (m) the predicted ejected mass during a nova 
outburst is ~ 10"^ Mq. 

It should be stressed that the amount of mass ejected during a nova outburst predicted 
by hydrodynamical models (~ 10^^ ^0) is systematically lower than the value observa- 
tionally determined in a small sample of nova systems, peaked around ~ 10^"^ Mq (see e.g. 
Delia Valle 2000 for a summary of results). However, the envelope masses inferred from 
observations are highly uncertain. In addition, the usual assumption that the ejected shells 
are almost homogeneously filled in (filling factor ~ 0.1-1) has been challenged by recent 
observations of the Nova T Pyx (Shara et al. 1997), suggesting values of the fiUing factor in 
the range 10~^-10"~^. Incomplete knowledge of the mass of nova ejecta is the most serious 
limitation to a quantitative evaluation of the role of novae as ^Li sources in the Galaxy. 

Crude estimates of the total amount of ^Li synthesized by Galactic novae have been given 
by Starrfield et al. (1978), Hernanz et al. (1996), and JH98. They all consistently show that 
novae can account for ~ 10% of the Galactic ^Li content. The contribution of novae to the 
Galactic evolution of ^Li has been considered by D'Antona & Matteucci (1991), Matteucci et 
al. (1995), Romano et al. (1999), in the framework of models of Galactic chemical evolution. 
However, a detailed study of the role of novae as ^Li producers necessarily suffers from several 
uncertainties. On one hand, the relevant quantity, i.e. the amount of ''Li ejected into the ISM 
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during a nova outburst, is the product of two poorly constrained factors: the ^Li abundance 
and the total mass ejected. On the other hand, the evaluation from first principles of the 
nova rate in the Galaxy requires the accurate knowledge of a number of physical processes 
and quantities which are neither theoretically nor observationally well determined, such as 
the white dwarf cooling timescale, the fraction of binary stars, and the fraction of binary 
stars that end up as a nova system. 

In this Section, we estimate the upper limit to the amount of ^Li synthesized by novae 
in the Galaxy, in particular deriving approximate analytical expressions for the abundance 
of ^Li as a function of the gas metallicity. In the framework of the "closed box" model for 
the solar neighborhood (see e.g. Tinsley 1980), it is easy to predict the evolution of the 
^Li mass fraction X-j[^) as a function of the gas fraction /i in the Galaxy. Let's define a nova 
mass ejection rate ip-^{t) = MejZ^N(^) (in Mq yr^"*^), where ^'N(t) is the nova rate (in yr^^), 
and assume that the rate of mass ejection by novae is proportional to the star formation 
rate, (pt^(t) = aip{t), with a independent of time. Indicating with {Xi)ej the average mass 
fraction of ''Li in nova ejecta, and Xrm the initial (cosmological) ^Li abundance, we obtain 

XM = + ^^[1 - /.«/(^-«)], (4) 

where R ~ 0.21 according to Galli et al. (1995) is the stellar returned fraction over the 
Galactic lifetime. The first term in eq. (H) represents the ^Li destruction by astration, the 
second term the ^Li production by novae. It is convenient to eliminate the gas fraction /x in 
favour of the metallicity of the gas Z, given by 

Pz 

ZifJ') = 'yZTr^^^' (5) 
where Pz — 7.9 x 10^^ (Galli et al. 1995) is the metal production factor. 

Notice that the contribution of novae to the evolution of ^Li in the Galaxy depends 
almost linearly, at late times, on the combination ^(Xy)^. The value of a, being constant, 
can be estimated at the present time t = tcab 

^ ^ Mej^N(tGal) ^ ^ ^ / ^ej \ f MtG.l) \ f V^(tGal) \ ^ 



Gal J 



where we have adopted the value of the nova rate recently proposed by Shafter (1997) and 
the star formation rate given by Metzger (1988). As for Mgj and {X'/)^] we have adopted the 
average values given by Romano et al. (1999) for t > 4.5 x 10'' yr, based on the results of 
JH98: Mej ^ 3 X 10"^ Mq and (X7)ej ^ 3 x 10"^ Clearly for a ~ 10"^ the role of novae in 
the Galactic evolution of ^Li is marginal, the ^Li abundance predicted by the simple model 
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at i = ^0 (^7 ~ 1 X lO-'^) being only ~ 10% of the meteoritic value (Xre ~ 9.3 x IQ-^). 
This result is in agreement with the order-of-magnitude estimates by Starrfield et al. (1978), 
Hernanz et al. (1996), and JH98. 

Several approximations have been made in the derivation of the results of this Section. 
However, they are hkely to lead to an overestimate of the predicted ^Li abundance. Consider 
for instance the assumption of Instantaneous Recycling Approximation (IRA) for novae. 
Owing to the long evolutionary timescales necessary to produce a nova system, the ratio a 
was considerably smaller in the past (see D'Antona & Matteucci 1991, Romano et al. 1999), 
and the approximation of a nova rate proportional to the SFR, with the proportionality 
constant a estimated at the present time, results in an overestimate of the past nova rate 
and therefore the ^Li abundance. Neglecting infall also results in an overestimate of the 
^Li abundance predicted by the model, if the infalling material has primordial composition 
and therefore acts as a diluting effect on the disk ^Li abundance. The same is true neglecting 
the delay of about ~ 1 Gyr estimated by D'Antona & Matteucci (1991) between the onset 
of star formation in the Galaxy and the birth of the first nova system. 

Finally, another independent argument against a dominant contribution of novae to 
the production of ^Li is the constraint on the abundance of isotopes hke ^^C, ^^N and ^^O 
produced copiously by novae according to JH98. The constraint that we discuss below is 
independent of any specific model of chemical evolution. In Table 2 we show the average 
mass fraction of these isotopes in nova ejecta, assuming the yields computed by JH98 for CO 
and ONe novae (third and fifth column), together with their abundances in the protosolar 
material (first column), representative of the ISM composition at t = (from Anders & 
Grevesse 1989, and Grevesse, Noels & Sauval 1996). In the fourth and sixth column we show 
the contribution of novae to the solar abundance of each element assuming that novae are 



TABLE 2 

Production of ^^C, ^^N, ^''O, and ''Li by novae 





^0 




Wej (a) 










3.6 X 10- 


-5 


1.0 X 10-^ 


100% (assumed) 


2.4 X 10-2 
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15N 


3.6 X 10- 


-6 


1.5 X 10-2 
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5.3 X 10-2 


220% 


1^0 


3. 1 X 10- 


-6 


1.1 X 10-2 


120% 


3.0 X 10-2 
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9.3 X 10" 


-9 


3.0 X 10-6 


12% 


9.2 X 10-^ 


14% 



(a) CO novae (JH98 models) 
(6) ONe novae (JH98 models) 
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the only producers of^^C. This is obviously not the case, since it is well known that low- and 
intermediate-mass stars produce significant amounts of this isotope (see e.g. Palla et al. 2000 
and references therein) . We see from Table 2 that even under this favourable assumption on 
the role of novae in the chemical enrichment of the Galaxy, their contribution to the solar 
Li abundance cannot be greater than ~ 10%. 

6. The production of Li by SNII 

Woosley et al. (1990) and Woosley & Weaver (1995, hereafter WW95) advanced the 
idea that production of ^Li via the so called neutrino process in SNII could account entirely 
for the Solar System ^Li abundance. Timmes, Woosley, & Weaver (1995), using a full grid 
of SNII models of various masses and mctallicities, predicted a lower ^Li production rate by 
the i/-process than Woosley et al. (1990), concluding that SNII contribute about one-half 
the solar '''Li abundance (see also Matteucci et al. 1995 on this point). 

With our model of GCE, assuming the ^Li yields of WW95, we obtain the results shown 
in Figure 7. Starting with an initial ^Li composition of Xyin = 1.0 x 10^^'^, wc found that 
SNII can account for ~ 40% of the meteoritic ^Li content, in agreement with Timmes et 
al. (1995), and Matteucci et al. (1995). In order to show the exact value of the metaUicity 
at which SNII in our GCE model start to contribute significantly, we also show in Fig. 7 the 
contribution of SNII computed with initial X-^i^ — 0. 

We should notice, however, that the SNII Li yields computed by WW95 have been 
questioned in recent theoretical studies of hydrodynamics and rotation in SNII (see e.g. 
Langer et al. 1999, and Heger, Langer, & Woosley 2000). In particular, Langer et al. (1999) 
showed that the inclusion of rotational mixing of the envelope of massive MS stars (later 
supposed to end as SNII) drastically reduces the amount of ^He present in the stellar interior. 
In non-rotating massive stars this isotope is found to be neither produced nor destroyed. As 
shown by Langer ct al. (1999) for a 15 Mq rotating model, ^Hc production factor is 0.1% 
of the correspondent case in a non-rotating WW95 model. Since ^Li is largely produced 
by the ^He(Q;,7)^Li reaction, the inclusion of rotation in massive star models is expected to 
drastically reduce the production of ^Li , with respect to the non-rotating cases. Thus, the 
results shown in Fig. 7 should be considered as upper limits for the contribution of SNII to 
the ^Li enrichment in the Galaxy. 
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7. The production of "Li by deep mixing in low-mass giants stars 

Pop. I low mass giant stars (M < 4 Mq) are thought to produce Li both in the RGB 
and AGB phases. More specifically, low mass stars, below 2.5 Mq, produce ^Li in the RGB 
and in the AGB by a deep mixing process such as explored by Wasserburg, Boothroyd, & 
Sackmann 1995, and Sackmann & Boothroyd 1999, hereafter SB99. This process is some- 
times called cool bottom processing. We prefer to refer to this as "deep mixing" because it 
is more indicative of the physics involved: the cool bottom of the envelope does not play any 
role. The nuclear processing occurs at the top of the H-shell and this is facilitated by mixing 
from the bottom of the convective envelope down to the deeper layers of the H-shell. This 
mechanism needs to be further explored to see if it can be efficient also in the mass range 
2.5 < M/Mq < 4.0. 

According to stellar evolution models, RGB stars should be characterized by a relatively 
low '^Li content: stars that leave the MS undepleted in ^Li , after the first-dredge up dihition 
has occurred, are expected to have ^Li abundances of the order of loge(^Li) ~ 1.5 — 1.9 
(Iben 1965, 1967a,b). Since most stars destroy ''Li on the MS, their ''Li content when they 
reach the RGB and in immediately subsequent phases should be much lower than the above 
values. Most field and cluster RGB stars show in fact very low ^Li abundances (—1.5 < 
loge(^Li) < 0), even lower than predicted by the models. An extra-mixing mechanism has 
been proposed to explain those low abundances (see e.g. Charbonnel, Brown, & Waller- 
stein 1998; Charbonnel & Balachandran 2000, and references therein). However, after the 
first discovery by Wallerstein &; Sneden (1982) of a Li-rich RGB (i.e. a RGB star with a 
''Li abundance higher than model predictions), several other Li-rich giants have been found 
(for recent surveys, see Castilho et al. 1998; Jasniewicz et al. 1999). Note that, although the 
number of presently known Li-rich RGB stars is relatively high, they represent only 1-2 % of 
all the giants with ^Li measurements (Wallerstein & Sneden 1982; Gratton & D'Antona 1989; 
Pilachowski, Hudek, & Sneden 1990; Pallavicini et al. 1990; Fekel & Balachandran 1993). 
Some of these Li-rich RGB stars have abundances even higher than the present ISM value 
(de la Reza & da Silva 1995; Balachandran 2000). 

Various suggestions have been made to explain the Li- rich giants with loge(^Li) > 2 (for 
stars with 1 < loge(^Li) < 2 a fresh '''Li production is not necessary since their ^Li abundance 
is consistent with the first dredge-up values). The high '''Li content of these giants may be 
related to external processes (Alexander 1967; Gratton & D'Antona 1989; Siess &: Livio 1999) 

or to internal processes such as the production of fresh ''Li (Fekel & Balachandran 1993; de la 
Reza, Drake, & da Silva 1996; SB99). As shown by SB99, ''Li can be created in low-mass RGB 
stars via the Cameron- Fowler mechanism associated with a deep-mixing below the convective 
envelope. This internal circulation (possibly driven by stellar rotation) transports envelope 
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material into the outer wing of the H-burning shell, where it undergoes partial nuclear 
processing, and then back to the envelope. SB99 choosed the free parameters of their model 
to match the low value of ^^C/^'^C observed in RGB stars. In order to produce the required 
additional ^^C, the advected material must reach temperatures high enough that ^He is 
burned, resulting in the creation of ^Be via ^He(a;,7)^Be. If extra-mixing is slow, '^Be is 
destroyed via ^Be(p,7)^B(e"'",i/)^Be or ^B(e~,z/)^Li, and any ^Li produced from ^Be electron 
captures is immediately burned up via ^Li(p,Q;)''He. However, for higher mixing speed ^Be 
can be transported out to cooler regions before electron capture takes place and the stellar 
envelope becomes enriched in ^Li . The production of ^Li in RGB stars is dependent on 
the mixing-speed, i.e. the stream mass flow rate. SB99 discussed the range of values of this 
mixing-speed and argued that it must be slower than the velocity of convection in RGB or 
AGB (~ 1 Mq yr^^), while the streams must move faster than the speed with which the 
H-shell burns its way outward. We show in Table 3 how the different values for the mixing- 
speed can influence the ^Li production (see also SB99 for more details on the 1 Mq model). 
It is important to notice that, as discussed by e.g. Charbonnel (1994) and Boothroyd & 
Sackmann (1999), the mass range in which the deep mixing is active is ~ 1.0-2.5 Mq. This 
is due to the fact that for low mass stars (< 2.5 Mq), the H-burning shell catches up to and 
erases the discontinuity while the star is still on the RGB; for higher masses (> 2.5 Mq) the 
star leaves the RGB before this can take place. 

We use the models by SB99 to estimate the contribution of low-mass RGB stars to the 
chemical evolution of the Galaxy. The ^Li production phase on the RGB is short (~ 10^ yr) 
compared to the total red giant lifetime (~ 5 x 10^ yr), and since the typical RGB mass-loss 
rate is rather low, they are not expected to contribute significantly to the ^Li enrichment 
in the ISM. On the other hand, during the TP- AGB phase, because of the higher mass-loss 
rate, the contribution to the ''Li ISM enrichment can be substantial. For these reasons, 
and since SB99 do not present spcciflc ^Li predictions for the AGB phase, in this work we 
make the following assumptions. We consider the SB99 ^Li predictions at the end of the 
RGB phase (when log(L/LQ) ~ 3.4) for 1 Mq at two different metallicities {Z — Zq and 
Z = 0.001) and for three different cases of mixing-speed (10-^ 10-^ 10"=^ Mq yr-^). To 
obtain the yield of ^Li from these models we multiply the ^Li abundances by the total mass 
ejected, e.g. for 1 Mq model we used Mej ~ 0.45 Mq (see e.g. Weidemann 1984 for the 
initial-flnal mass relationship). We assume a mass range of 1-2.5 Mq, in agreement with 
the above discussion, and since the authors showed only the results for the 1 Mq model, we 
interpolate the ^Li production by deep mixing in the mass range 1 Mq - 2.5 Mq. To derive 
the 2.5 Mq ^Li yields we follow SB99 assuming that all ^He is converted in ^Li , and that 
the abundance of ^He in stars scales as M~^ (Schatzman 1987). The ^Li yields obtained 
for 1 Mq and 2.5 Mq, with different metallicities and different mixing speed, are shown in 
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Table 3. In our calculations, in order to derive an upper limit to the Li production, we 
adopt the highest ^Li abundance predicted by the SB99 models, i.e. that obtained with the 
model with 10~^ Mq yr~^ mixing speed. 

The contribution of low-mass giant stars to the Galactic evolution of ^Li is shown in 
Figure 8. These stars enrich ^Li up to Xj ~ 2.2 x 10~^ at the epoch of formation of the Sun, 
i.e. ~ 24% of the solar value. Again, another more efficient source is required to explain the 
ISM ^Li abundance. As in the case of SNII (Fig. 7) wc show for comparison the contribution 
from these stars to ''Li enrichment starting with a zero initial ''Li composition. 

8. Conclusions 

We have addressed the long-standing problem of the Galactic chemical evolution of 
^Li with the aim of clarifying the role of the different stellar contributions. Wc have ana- 
lyzed four possible stellar sources of ''Li : SNII, novae, low-mass giants and IMS-AGB stars. 
For low-mass giants, novae, and SNII we have critically examined the available ^Li yields 
and discussed the possible extrapolations when no model predictions were available. In the 
case of SNII and novae we re-examined the ''Li yields in the light of recent hydrodynamical 
computations. For IMS-AGB stars we presented and used here new results for nucleosyn- 
thesis calculations based on evolutionary AGB models by Frost (1997). In particular, we 
discussed how the interplay between the HBB process and a phase of high mass-loss before 
the evolution off the AGB may constitute a key process for ''Li enrichment in the Galaxy. 
Although we are not yet able to quantify the contribution of IMS-AGB stars, we have ex- 
plored different realistic possibilities in terms of mass-loss rate and mass range for the HBB 
process. 

The work presented here is summarized in Figure 9, where we compare the results of 
our model of chemical evolution with the sample of observational data presented in Fig. 1. 
The predicted ^Li abundance resulting from all stellar sources considered in this paper is 
plotted in Fig. 9 vs. metallicity (upper panel) and vs. time (lower panel). For IMS-AGB we 
used a VW93 mass-loss increased by a factor of 50 (as discussed in § 4.2), and we have also 
taken into account the results obtained for the 3 Mq AGB model presented in § 4.3. 

As we discussed in Sect. 2.1, we also ran our model with an initial ''Li abundance higher 
by a factor 2 than the Spite plateau value. The result is shown for comparison in Fig. 9. 
Notice that the predicted ''Li abundance at the time of formation of the Sun is virtually 
the same in both cases. In the lower panel of Fig. 9 we plot the individual contributions 
from the different stellar sources, with the exception of the contribution from SNII since 
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we believe that the inclusion of rotation and mixing in the supernova models leads to a 

drastic reduction of the ^Li yields (see § 6). Note that, from an observational point of 
view, it is critical to determine the metallicity where the rise of the ^Li abundance from the 
Spite plateau occurs; the available data indicate that the rise from the plateau occurs at 
metallicities [Fe/H] between ~ —1 and —0.8, but additional data in the metallicity range 
— 1 < [Fe/H] < —0.3 are needed in order to better constrain this value. 

Our main conclusions are the following: 

(i) In the light of the available stellar models, small contributions to the meteoritic 
^Li abundance come from novae (~ 10% or less, see § 5) and low-mass giant stars (~ 20%, 
see § 7). As for SNII, we predict a contribution to the solar system ^Li abundance of ~ 40% 
with the standard yields by Woosley & Weaver (1995), and a contribution less than 10% if 
the results of the latest hydrodynamical simulations of the supernova explosion are taken 
into account (see § 6). Spallation reactions between ISM and cosmic-ray nuclei can provide 
an additional at least ~ 10-20% of the ^Li abundance in the solar system (see e.g. Lemoine 
et al. 1998). 

(ii) Figure 9 shows that even with our best-fit model we are not able to reproduce the 
meteoritic abundance, although both cluster and ISM abundances are fitted fairly well. Our 
best-fit model is in good agreement with the observed loge(^Li) vs. [Fe/H] distribution 
up to [Fe/H] ~ —0.4. This indicates that the metallicity at which the rise from the Spite 
plateau occurs is consistent with the [Fe/H] values at which the IMS-AGB stars contribute 
to the chemical enrichment of the ISM. We remind that the observational data for [Fe/H] 
> —0.4 are mostly taken from Balachandran (1990) (see discussion in § 2); nevertheless, 
as the comparison between photometric and spectroscopic metallicities shows, part of these 
stars might have higher [Fe/H] values. 

(Hi) Our best-fit model and, more specifically, the contribution from IMS-AGB, is based 
on two major assumptions: first the high mass-loss phase at the end of the evolution of 
these stars must start earlier on the AGB than the standard predictions by VW93. This is 
supported by recent IR observations, discussed in § 4. Second, we also need a contribution 
to ^Li through HBB production from stars with 3 Mq < M < 4 Mq and low metallicities 
(see preliminary model results for M < 4 Mq in Lattanzio et al. 2001). Should these two 
assumptions be not valid, the contribution from intermediate-mass AGB would be much 
lower. We stress however that, since neither novae, nor SNII, nor low-mass giant stars 
seem to produce enough ^Li to account for the present-day abundance, we believe that 
intermediate-mass AGB stars remain at present the best candidates as ^Li factories. 

Our results and conclusions are strongly based on the AGB models discussed in § 4 
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(see also Frost 1997). A different view has been recently expressed by Ventura et al. (2000). 
They presented a grid of IMS-AGB models of different masses and tested the sensitivity of 
^Li yields to different mass-loss rate prescriptions on the basis of a comparison with AGB 
stars in the Magellanic Clouds. Ventura et al. (2000) concluded that, with their mass-loss 
calibration, IMS-AGB stars do not contribute significantly to the ^Li enrichment of the ISM. 
As they properly noticed, details of ''Li production depend on the input parameters of the 
stellar model, mainly the treatment of convection. In addition, Ventura et al. (2000) focused 
their work on the analysis of the strength of the mass-loss rate during the AGB phase. 
We instead followed a different approach, testing different times during the AGB phase at 
which the high mass-loss rate starts, and looking for the consequences on ^Li yields. As 
we demonstrated in § 4, anticipating by few thermal pulses the beginning of the superwind 
phase can have significant impact in the mass of ^Li ejected. 

Work in progress includes an extension of the grid of AGB models to lower masses 
(M < 4 Mq) and lower metallicities {Z < 0.004) (see prehminary results in Lattanzio et 
al. 2001), in order to analyze the possibility of ^Li production via HBB in these stars. Another 
interesting point that we just mentioned in this paper, is a re- analysis of the contribution 
of GCR to the production of ^Li (Valle et al., in preparation). Finally, new ^Li data for a 
statistically significant sample of stars in the critical metallicity range —1 < [Fe/H] < —0.3 
are being analyzed; the data will provide stringent observational constraints on the Galactic 
^Li abundance in a metallicity range where IMS-AGB stars are expected, on the basis of this 
work, to give their main contribution. 

We thank S. Shore for interesting discussions and comments on the manuscript. C.T. 
also thanks R. Gallino for useful suggestions and encouragements during the development of 
this work. We also thanks the anonymous referee for very usefuU comments that improved 
our paper. The research of C.T., S.R. and D.G. is partially supported by grants COFIN98 
and COFIN2000. 
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TABLE 3 

'''Li yields (in solar masses) by deep mixing in low-mass giant stars 



M 


mixing speed 




Z = 0.001 


(Me) 


(Mo yr-i) 






1.0 


10-3 


3.8 X 10-1° 


5.1 X 10-1° 




10-^ 


2.0 X 10-9 


7.6 X 10-1° 




10-^ 


1.0 X 10-^ 


3.6 X 10-1° 


2.5 


10-3 


2.1 X 10-^° 


3.0 X 10-1° 




10-^ 


7.6 X 10-^° 


6.7 X 10-1° 




10-'' 


6.0 X 10-" 


2.1 X lO-'^ 



(a) SB99 models. 

(6) Our extrapolations from (a). 
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Fig. 1. — Upper panel: loge(^Li) vs. [Fe/H] for field stars. Observations are from Full- 
bright (2000) {filled circles), Bonifacio & Molaro (1997) {open squares), Balachandran (1990) 
{open triangles), Lambert et al. (1991) {open circles). The box represents the region occupied 
by stars in open clusters with undepleted ^Li . The meteoritic and the (solar) photospheric 
values are also shown by a dotted circle and an asterisk, respectively. Lower panel: logeC^Li) 
vs. age for Galactic open clusters {filled triangles), and globular clusters {box). We also 
include different measurements of the ISM value {filled circles), the meteoritic value {dotted 
circle), and the cosmological ^Li abundance adopted in this work (see text for references). 
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Fig. 2. — Surface ^Li abundance in the envelope of 4. 5. 6 Mq stars with Z = 0.004, 0.008, 
0.02, vs. time on the AGB. Each circle represents a thermal pulse in the model. 
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g. 3. — Surface C/0 ratios vs. ^Li surface abundances for the nine AGB models described 
the text. Symbols are as Fig. 2. 
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Fig. 4. — Surface ^Li abundance vs mass-loss rate (in Mq yr~^) for the nine AGB models 
described in the text, following the mass-loss prescriptions of Vassiliadis & Wood (1993). 
Symbols are as Fig. 2. 
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Fig. 5. — Surface ^Li abundance vs. stellar mass for the AGB models discussed in the text. 
The open circles are for the Vassiliadis & Wood (1993) mass-loss, the open triangles are for 
the Vassiliadis & Wood (1993) mass-loss increased by a factor 50. Each open circle and open 
triangle represents a thermal pulse. 
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Fig. 6. — Contribution of IMS-AGB stars to the Galactic evolution of ^Li , for two different 
choices of the IMS-AGB mass range: 4-6 Mq (upper panel), and 4-6 Mq for Z > 0.004 and 
3-6 Mq with Z < 0.004 (lower panel) (see also text for details). In both cases wc adopt 
the mass-loss rate prescription of Vassiliadis & Wood (1993) with modification (sec text). 
Symbols are as in Fig.l (upper panel). Lines refer to the GCE model results for halo (dotted), 
thick disk (dashed), and thin disk (solid). Thick lines show the ^Li evolution obtained with a 
zero initial abundance of ^Li , to emphasize the contribution of IMS-AGB stars in the thick 
disk (dashed line) and in the thin disk (solid line). 
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Fig. 7. — Contribution of SNII to the Galactic evolution of ^Li (line types are as in Fig. 6). 
Thick lines show the ^Li evolution obtained with a zero initial abundance of ^Li , to emphasize 
the contribution of SNII in the thick disk [dashed line) and in the thin disk [solid line). 
Symbols are as in Fig.l (upper panel). 
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Fig. 8. — Contribution of low-mass giant stars via deep mixing process to the Galactic 
evolution of ^Li (line types are as in Fig. 6). Thick lines show the ^Li evolution obtained 
with a zero initial abundance of ^Li , to emphasize the contribution of low-mass giants in 
the thick disk {dashed line) and in the thin disk {solid line) Symbols are as in Fig.l (upper 
panel) . 
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Fig. 9. — Upper panel: evolution of ^Li as a function of [Fe/H] according to our model taking 
into account the contribution of novae, low-mass giants and IMS-AGB stars (line types are 
as in Fig. 6). We also show the results obtained starting with an initial ^Li abundance 
higher by a factor 2 than the Spite plateau. Lower panel: log e(^Li) vs. age for Galactic open 
clusters {filled triangles) and globular clusters (same symbols as in Fig. 1). We indicate the 
individual contributions of novae {long-dashed line), low-mass giants {short-dashed line), and 
IMS-AGB stars {dotted line). The solid line shows the total contribution from all sources 
considered. 



